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Abstract 

The gluon polarisation in the nucleon has been determined by detecting charm 
production via meson decay to charged K and tt in polarised muon scattering off 
a longitudinally polarised deuteron target. The data were taken by the COMPASS 
Collaboration at CERN between 2002 and 2006 and corresponds to an integrated 
luminosity of 2.8 fb~^. The dominant underlying process of charm production is 
the photon-gluon fusion to a cc pair. A leading order QCD approach gives an 
average gluon polarisation of {Ag/g)^= —0.49 ± 0.27(stat) ± O.ll(syst) at a scale 
/i^ ~ 13(GeV/c)^ and at an average gluon momentum fraction (x) ~ 0.11. The 
longitudinal cross-section asymmetry for production is presented in bins of the 
transverse momentum and the energy of the meson. 



PACS: 13.60.-r, 13.88.+e, 14.20.Dh, 14.70.Dj 

keywords: Inelastic muon scattering; Spin; Asymmetry; Gluon polarisation 



submitted to Physics Letters 



The COMPASS Collaboration 



M. Aleksccv'^^) , V.Yu. Alcxakhin^), Yu. Alexandrov^'^) , G.D. Alexeev^), A. Amoroso^*), 
A. Austrcgisilioii'18), B. Baddek^i), F. Balestra^s), J. BalF), J. Barth^), G. Baum^), 
Y. Bedfer23), J. Bernhard^^), R. Bertini^^), M. Bettinelli^^) , R. Birsa^^^, J. BisplinghofF^) , 
P. Bordalo^^'^\ F. Bradamante^^) , A. Bravar^^), A. Bressan^^), Q. Brona^^\ E. Burtin^^), 
M.P. Bussa^s), A. Chapiro^^), M. Chiosso^^), S.U. Chung^^), A. Cicuttin25'27)^ m. Colantoni^^) , 
M.L. Crespo25'27)^ s. Dalla Torre^^), T. Dafni^^), S. Das^), S.S. Dasgupta^), O.Yu. Dcnisov^^'b) ^ 

L. Dhara'^), V. Diaz25'27)^ A.M. Dinkclbachi^) , S.V. Donskov^^), N. Doshita^''^^) , V. Duic^^), 
W. Diinnweber^''), A. Efremov®), A. El Alaoui^^^P.D. Eversheim^) , W. Eyrich^\ M. Faessler^^\ 
A. Ferrero^^'^^), M. Finger^"), M. Finger jr.^), H. Fischer^^), C. Franco^^), J.M. Friedrich^^\ 
R. Garfagnini^s), F. Gautheron^), O.P. Gavrichtchouk^) , R. Gazda^^), S. Gerassimov^'^'^^), 
R. Geyeri^), M. Giorgi^e), B. Gobbo^s), S. Goertz^.^), S. Grabmiilleris) , O.A. Grajek^i), 
A. Grasso^s), B. Grubc^^), R. Gushtcrski*^) , A. Guskov^), F. Haas^^), R. Hagemann^°),D. von 
Harrach^"^), T. Hasegawa^^-* , J. Heckmann^), F.H. Heinsius^'^'*, R. Hermann^^\ 

F. Herrmann^°),C. Hefi^^, F. Hinterberger^\ M. von Hodenberg^^^N. Horikawa^^'^\ 
Ch. H6ppneri8),N. d'Hose^^), C. Ilgner^^'^'^) , S. Ishimotoi9'<^) , O. Ivanov^), Yu. Ivanshin^), 

T. Iwata^^), R. Jahn^), P. Jasinski^^), G. Jegou^^), R. Joosten^), E. Kabufi^^), 
W. Kaferio),D. Kang^"), B. Ketzer^^), G.V. Khaustov^^), Yu.A. KhokhW^), 
J. Kieferio),Yu. Kisselev^-^), F. Klein*^), K. Klimaszewski^i) , S. Koblitz^''), J.H. Koivuniemi^) , 
V.N. Kolosov^^) , E.V. Komissarov^'+) , K. Kondo^'^^^ , K. K6nigsmann^°) , I. Konorov^^'^^^ , 

V.F. Konstantinov^^) , A. Korzenev-*^^'^^ , A.M. Kotzinian^'^^),0. Kouznetsov^'^^^ , 
K. Kowalik^^'22\ M. Kramer^^), A. KraP^), Z.V. Kroumchtein^) , R. Kuhn^^), F. Kunne^^), 
K. Kurek^i), J.M. Le GofP^), A.A. Lednev22), A. Lehmann^), S. Levorato^e), 
J. Lichtenstadt^^) , T. Liska^^), A. Maggiora^^), M. Maggiora^^) ,A. Magnon^^), G.K. Mallot^^), 
A. Mann^^), C. Marchand^^), J. Marroncle^^) , A. Martin^^), J. Marzec^^), F. Massmann^), 
T. Matsuda^^) , A.N. Maximov^'+) ,W. Meyer^) , T. Michigami^^) , Yu.V. Mikhailov^^) , 
M.A. Moinestcr24) , A. Mutter lO'^"^), A. Nagaytscv^), T. Nagel^^), J. Nassalski^^) , S. Negrini^), 

F. Ncrlingio), S. Neubert^^), D. Neyret^^), V.I. Nikolaenko^^) , A.G. Olshcvsky^), 
M. Ostrick^'^"^), A. Padee'^^\ R. Panknin*^^ S. Panebianco^^),D. Panzieri^°\ B. Parsamyan^^\ 
S. Paul^^\ B. Pawlukiewicz-Kaminska^^\ E. Perevalova^\ G. Pesaro^^^ D.V. Peshekhonov^^ , 

G. Piragino^s), S. Platchkov^^) , J. Pochodzalla^^) , J. Polaki2,26)^ y a. Polyakov^^), 

G. Pontecorvo*), J. Pretz^), C. Quintans^^), J.-F. Rajottc^^), S. Ramos^-'^'^) , V. Rapatsky^), 

G. Reicherz^), D. Reggiani"), A. Richter'^), F. Robinct^''^) , E. Rocco^^), E. Rondio^^), 
D.I. Ryabchikov^^) , V.D. Samoylenko^^) , A. Sandacz^^), H. Santos^^-'') , M.G. Sapozhnikov^^ , 
S. Sarkar^), LA. Savin^), G. Sbrizza^^), P. Schiavon^^), C. Schill^o), L. Schmitt^^'^) , 
W. Schroder^), O.Yu. Shevchenko^) , H.-W. Siebert^^), L. Silva^^), L. Sinha^), A.N. Sissakian^), 
M. Slunecka^), G.I. Smirnov^), S. Sosio^s), F. Sozzi^^), A. Srnka^), M. Stolarski^i'ii), 
M. Sulci2), R. Sulej32), S. Takekawa^e), S. Tcssaro^^), F. Tcssarotto^^) , A. Tcufel^), 
L.G. Tkatchev^), G. Venugopal^), M. Virius^^), N.V. Vlassov^), A. Vossen^o), Q. WeitzeF), 
K. Wenzli°),R. Windmolders^) , W. Wislicki^^^) , H. Wollny^o), K. Zaremba^^), M. Zavertyaev^^) , 
E. Zemlyanichkina^) , M. Ziembicki^^) ^ j Zhao-*^^'^^) , N. Zhuravlev^) and A. Zvyagin^''') 



Universitat Bielefeld, Fakultat fiir Physik, 33501 Bielefeld, Germany^) 

Universitat Bochum, Institut fiir Experimcntalpliysik, 44780 Bochum, Germany^) 
Universitat Bonn, Helmholtz-Institut fiir Strahlen- und Kernphysik, 53115 Bonn, Germany^^ 
Universitat Bonn, Physikalisches Institut, 53115 Bonn, Germany^) 
Institute of Scientific Instruments, AS CR, 61264 Brno, Czech Republic^^ 
Burdwan University, Burdwan 713104, India*^-* 

Matrivani Institute of Experimental Research & Education, Calcutta-700 030, India'-* 
Joint Institute for Nuclear Research, 141980 Dubna, Moscow region, Russia 
Universitat Erlangen-Niirnberg, Physikalisches Institut, 91054 Erlangen, Germany*' 
Universitat Freiburg, Physikalisches Institut, 79104 Freiburg, Germany*' 
CERN, 1211 Geneva 23, Switzerland 

Technical University in Libcrcc, 46117 Liberec, Czech Republic^' 
LIP, 1000-149 Lisbon, Portugal^) 

Universitat Mainz, Institut fiir Kernphysik, 55099 Mainz, Germany*' 

University of Miyazaki, Miyazaki 889-2192, Japan*^' 
Lebedev Physical Institute, 119991 Moscow, Russia 

Ludwig-Maximilians-Universitat Miinchen, Department fiir Physik, 80799 Munich, Germany*'^' 
Technische Universitat Miinchen, Physik Department, 85748 Garching, Germany*''' 
Nagoya University, 464 Nagoya, Japan"^' 

Charles University, Faculty of Mathematics and Physics, 18000 Prague, Czech Republic^' 
Czech Technical University in Prague, 16636 Prague, Czech Republic^' 

State Research Center of the Russian Federation, Institute for High Energy Physics, 142281 Protvino, 

Russia 

CEA DAPNIA/SPhN Saclay, 91191 Gif-sur-Yvette, France 

Tel Aviv University, School of Physics and Astronomy, 69978 Tel Aviv, Israel™' 

Trieste Section of INFN, 34127 Trieste, Italy 

University of Trieste, Department of Physics and Trieste Section of INFN, 34127 Trieste, Italy 

Abdus Salam ICTP and Trieste Section of INFN, 34127 Trieste, Italy 

University of Turin, Department of Physics and Torino Section of INFN, 10125 Turin, Italy 
Torino Section of INFN, 10125 Turin, Italy 

University of Eastern Piedmont, 1500 Alessandria, and Torino Section of INFN, 10125 Turin, Italy 
Soltan Institute for Nuclear Studies and University of Warsaw, 00-681 Warsaw, Poland"' 
Warsaw University of Technology, Institute of Radioelectronics, 00-665 Warsaw, Poland"' 
Yamagata University, Yamagata, 992-8510 Japan*^' 

Deceased 

Also at 1ST, Universidade Tecnica de Lisboa, Lisbon, Portugal 

On leave of absence from JINR Dubna 

Also at Chubu University, Kasugai, Aichi, 487-8501 Japan^' 

Also at KEK, 1-1 Oho, Tsukuba, Ibaraki, 305-0801 Japan 

Also at GSI mbH, Planckstr. 1, D-64291 Darmstadt, Germany 

Supported by the German Bundcsministcrium fiir Bildung und Forschung 

Suppported by Czech Republic MEYS grants ME492 and LA242 

Supported by DST-FIST II grants. Govt, of India 

Supported by the Shailabala Biswas Education Trust 

Supported by the Portuguese FCT - Fundacao para a Ciencia e Tecnologia grants 
POCTI/FNU/49501/2002 and POCTI/FNU/50192/2003 

Supported by the MEXT and the JSPS under the Grants No. 18002006, No.20540299 and No. 18540281; 
Daiko Foundation and Yamada Foundation 

Supported by the DFG cluster of excellence 'Origin and Structure of the Universe' ( www. uni verse- 
cluster, de) 

Supported by the Israel Science Foundation, founded by the Israel Academy of Sciences and Human- 
ities 

Supported by Ministry of Science and Higher Education grant 41/N-CERN/2007/0 
Supported by KBN grant nr 134/E-365/SPUB-M/CERN/P-03/DZ299/2000 



1 



1 Introduction 

Pioneering experiments on the spin structure of the nucleon performed in the sev- 
enties at SLAC [1] were followed by the EMC experiment at CERN which obtained a 
surprisingly small quark contribution to the proton spin [2], in contrast to the naive ex- 
pectation that the spin of the nucleon is built mainly from valence quark spins [3]. This 
result triggered extensive studies of the spin structure of the nucleon in polarised lepton 
nucleon scattering experiments at CERN by the SMC [1] and COMPASS [5], at SLAC 
[6], at DESY [7] and at JLAB [8] as well as in polarised proton-proton collisions at RHIC 
P, [in]. As a result, the parton helicity distributions in the nucleon were extracted using 
perturbative QCD analyses. The contribution of the quark spins to the nucleon spin is 
now confirmed to be around 30%, smaller than 60%, the value expected from the EUis- 
Jaffe sum rule |Tl]. The reduction from the naive expectation of 100% can be explained 
by the relativistic nature of quarks (e.g. in the MIT bag model) [12]. However, due to the 
limited range in the four-momentum transfer squared, Q^, covered by the experiments, 
the QCD analyses (e.g. [5]) show limited sensitivity to the gluon helicity distribution 
as a function of the gluon momentum fraction x, Ag{x), and to its first moment, AG. 
(The perturbative scale, /i^, in these QCD analyses is set to Q^.) The determination of 
Ag{x) from QCD evolution has therefore to be complemented by direct measurements in 
dedicated experiments. 

The average gluon polarisation in a limited range of x, (Ag/g)^, has been determined 
in a mo del- dependent way from the photon-gluon fusion (PGF) process by HERMES [13j, 
SMC [13] and COMPASS [15]. These analyses used events containing hadron pairs with 
high transverse momenta, px, (typically 1 to 2 GeV/c) with respect to the virtual photon 
direction. PYTHIA [16j was used by HERMES and by COMPASS for the analysis of small 
events, while LEPTO [llj was used in SMC and the ongoing COMPASS analysis for 
> 1 (GeV/c)^ events. This method provides good statistical precision but relies on 
Monte Carlo generators simulating QCD processes. The measurements point towards a 
small value of the gluon polarisation at x ~ 0.1. This is in line with recent results from 
PHENIX [9j and STAR [ID] at RHIC. 

Taking into account quark and gluon orbital angular momenta, L, the nucleon spin 
projection (in units of h) is 

5, = ^ = ^AE + AG + L,, (1) 

where AE is the first moment of the sum of the quark helicity distributions. The decom- 
position of Eq. (Q, however gauge dependent, is defined in the infinite momentum frame 
where the quark parton model is valid. 

Here we present a new result on (Ag/g)^ from muon-deuteron scattering 1^ The 
gluon polarisation is determined assuming that open-charm production is dominated by 
the PGF mechanism yielding a cc pair which fragments mainly into D mesons. This 
assumption is supported by the measurements of F2 in the COMPASS kinematic domain 
[18] and further discussed in [19]. This method has the advantage that in lowest order of 
the strong coupling constant there are no other contributions to the cross-section; however, 
it is statistically limited as will be shown in section [3l In the present analysis only one 
charmed meson is required in every event. This meson is selected through its decay in one 

The present result includes a larger data sample and an improved analysis method and thus supersedes 
the one given in Ref. [20] , 
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of the two channels: D*(2010)+ ^ D'^n^^^^ R-n+n^^^^ {D* sample) and D° K-n+ 
{D^ sample) and their charge conjugates. 

2 Experimental set-up 

The data were collected between 2002 and 2006 with the COMPASS experiment at 
the M2 muon beam line of the CERN SPS. A detailed description of the experiment for 
the years 2002 to 2004 can be found in Ref. [21]. For the 2006 data taking the polarised 
target and the spectrometer were considerably upgraded. 

The measurements were performed using a /i"*" beam of 160 GeV/c. The beam muons 
originating from vr"*" and decays are naturally polarised with an average polarisation, 
P^, of about 80% with a relative uncertainty of 5% [22]. The momentum of each incoming 
muon is measured upstream of the experimental area with a precision of Ap/p < 1% in a 
beam momentum station consisting of layers of scintillators. The incoming muon direction 
and position is measured with a detector telescope in front of the target. A precision of 
30 /irad is obtained for the track direction. 

The polarised ^LiD target is housed in a superconducting solenoid with a polar angle 
aperture of 70 mrad in 2002 to 2004. The target consisted of two 60 cm long cells (upstream 
u, downstream d), separated by 10 cm, longitudinally polarised with opposite orientations. 
The spin directions were reversed every eight hours by rotating the field of the target 
magnet system. The target was upgraded in 2006 with a new solenoid with an aperture of 
180 mrad. To reduce the systematic errors due to the different spectrometer acceptances 
for the upstream and downstream cells, a 3-cell target configuration was installed. A 
central 60 cm long cell is placed in-between two 30 cm long cells with polarisations opposite 
to the central oneEl In this set-up the average acceptances for both spin directions are 
very similar and therefore the magnetic field was rotated only once per day. The average 
target polarisations, Pt, were 50% with a relative uncertainty of 5%. The dilution factor 
/, accounting for the fraction of polarisable nucleons in the target, is about 0.4, since 
the ®Li nucleus basically consists of a ^He core plus a deuteron. The exact value of / is 
kinematics dependent and is calculated as described in Ref. p3]. Its relative uncertainty 
is 5%. 

The two-stage COMPASS spectrometer is designed to reconstruct the scattered 
muons and the produced hadrons in a wide momentum range. Particle tracking is per- 
formed using several stations of scintillating fibres, micromesh gaseous chambers and 
gas electron multiplier chambers for the small angles tracks. Large area tracking devices 
comprise gaseous detectors (drift chambers, straw tubes and multiwire proportional cham- 
bers). The detectors are placed around the two spectrometer magnets. The direction of 
the tracks reconstructed at an interaction point in the target is determined with a preci- 
sion better than 0.2 mrad and the momentum resolution for charged tracks detected in 
the first spectrometer is about 1.2% whereas is it about 0.5% in the second spectrometer. 
The achieved longitudinal vertex resolution varying from 5 mm to 25 mm along the target 
allows assigning each event to a particular target cell, i.e. a specific target spin direction. 
For 2006 the tracking systems in the first stage were adapted to match the increased 
aperture of the polarised target magnet. The trigger is formed by several hodoscope sys- 
tems supplemented by two hadron calorimeters. Muons are identified downstream of the 
hadron absorbers. A Ring Imaging CHerenkov counter (RICH) with a C4F10 radiator is 
used in the first spectrometer stage for charged particle identification. It is equipped with 

■^^ In 2006 u and d stand for the central target cell and for the sum of the outer target cells, respectively. 
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multiwire proportional chambers with Csl photocathodes to detect the UV Cherenkov 
photons. The RICH, too, underwent a considerable upgrade for the 2006 data taking. 
In the central part, the photon detectors were replaced by multi-anode photomultiplier 
tubes, yielding considerably higher photon detection efficiency along with a much faster 
response. For the outer parts the readout electronics was refurbished, allowing a signifi- 
cant reduction of the background. The data taking amounted to 40 weeks in 2002 to 2006 
and corresponds to an integrated luminosity of 2.8 fb^^. 

3 Data selection 

In the present analysis the selection procedure required an incoming muon, a scat- 
tered muon, an interaction vertex in the target and at least two additional tracks. The 
kinematic variables like the four-momentum transfer squared Q^, the relative energy trans- 
fer and the Bjorken variable xbj = Q'^ /2MEy, where M is the nucleon mass and E 
the incident muon energy, are calculated from the four-momenta of the incident and 
scattered muon. No kinematic cuts are applied on Q"^, y or xbj. Thus the selected data 
sample includes the events with an interaction vertex from quasi-real photo-production 
Q"^ ^ m^?/^/(l — to a Q"^ of about 100 (GeV/c)^. Note that all the events are in the 
deep inelastic region, i.e. the invariant mass of the final state, W , is larger than 4 GeV/c^. 

The mesons are reconstructed through their Kir decay which has a branching 
ratio of 3.9%. Due to multiple Coulomb scattering of the charged particles in the solid 
state target the spatial resolution of the vertex reconstruction is not sufficient to separate 
the production and decay vertices. The mesons are selected using the invariant 
mass of their decay products. 

To reduce the large combinatorial background only identified Kt[ pairs are used. The 
identification in the RICH starts from reconstructed tracks with measured momenta. The 
likelihood for different mass hypotheses and for a background hypothesis are computed 
for each track, using the angles between the track and the detected Cherenkov photons. 
The likelihood functions, used in this computation, were defined from the corresponding 
expected angular distribution of photons; the expected distribution for background was 
obtained using a sample of photons not associated to reconstructed tracks. Particles are 
identified as kaons or pions on the basis of the likelihood associated to the pion, kaon, 
proton and background hypotheses. The procedure restricts the studied events to a sample 
with at least one kaon and one pion of momenta exceeding the Cherenkov threshold of 

9.1 GeV/c and 2.5 GeV/c, respectively. Simulations using the AROMA [27J generator and 
a full spectrometer simulation based on GEANT have shown that about 70 % (90%) of 
kaons (pions) coming from decays exceed this threshold for the reconstructed sample. 

All events have to satisfy a kinematic cut: z > 0.2, where z is the fraction of the 
energy of the virtual photon carried by the meson candidate. They are further divided 
into a D* and a sample, analysed independently. In the former one an additional track 
with a proper charge, a slow pion candidate, is demanded at the vertex. RICH likelihoods, 
used to reject electrons from those candidates, reduce the combinatorial background by a 
factor two. Furthermore, in the case of the D*, a cut on the mass difference is imposed, 

3.2 MeV/c2 < Mk^^^^^^ - Mk^ - < 8.9 MeV/c^, where Mx..,,_ and M^, are the 
masses of the D* and the candidates, respectively. Finally it was demanded that 
|cos6'*| < 0.9 for the D* sample and |cos^*| < 0.65 for the D^, where 6* is the decay 
angle in the D'^ cm. system relative to the flight direction. The events entering the 
D* sample are not used in the sample. The resulting mass spectra for the D'^ and D* 
samples with one Ktt pair in the mass range —400 MeV/c^ < M^n — M^o < 400 MeV/c^ 
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are displayed in Fig. [H A signal to background ratio in the signal region of about 1 is 
obtained for the D* sample and of about 0.1 for the sample with a mass resolution of 
about 22 MeV/c^ and 25 MeV/c^, respectively. The number of mesons is about 8,700 
and 37,400 in the D* and the samples. 

For the final event samples the mean value of is 0.65 (GeV/c)^, xbj ranges from 
1 ■ 10^^ to 0.6 with a mean value of 0.04 and y from 0.1 to 1 with a mean value of 0.55. 
Note that the perturbative scale for the selected events is not given by Q^, but by the 
transverse mass of the charmed quarks, = 4(m^ +Pt)- 

4 Method 

This section describes the determination of the gluon polarisation from the event 
samples collected in two different spin configurations and target cells. The same method 
is used in section [6] for the asymmetry determination. The number of events collected in 
a given target cell and time interval is 

l + wf^A^^^'^'^^^ + ^AB')! . (2) 
\s + b s + b J 

Here, ^m^^m'^^x _ _ o-TT^^'j^Q-Ti _|_ o-TT^^ where the arrows indicate the relative beam 
and target spin orientations, is the longitudinal double spin cross-section asymmetry of 
the events in the central peak of Fig. [T]and A-q is the corresponding asymmetry originating 
from the combinatorial background events in the mass spectra. Also, m = M^n, and X 
denote a set of kinematic variables describing an event (Q^, y, z...), while a, (p and n are 
the spectrometer acceptance, the integrated incident muon flux and the number of target 
nucleons, respectively. The differential unpolarised cross-sections of signal and background 
events folded with the experimental resolution as a function of m and X are represented 
by s = s{m,X) and b = b{m,X), respectively. The ratio s/{s + b) will be called "signal 
purity". In the present analysis the background is a combinatorial background and the 
signal purity can be extracted from the data using the invariant mass distributions of 
Fig. [H This is in contrast to the high-px analyses, where the physical background has to 
be estimated using a Monte Carlo simulation (MC) [131 [IH HH] ■ Information on the gluon 
polarisation is contained in A^^^^ ^ ^ which can be decomposed in LO QCD as 

A'^^^'^'^°^(X) = aLL(X)^(X). (3) 

9 

Here oll is the analysing power of the fig — > fi'cc process which includes the so-called 
depolarisation factor D accounting for the polarisation transfer from the lepton to the 
virtual photon. The background asymmetry Ab can be written as the product of the 
virtual photon asymmetry and the depolarisation factor = DA^ and is assumed to 
be independent of m. 

In the present analysis the average gluon polarisation {Ag/g),^ and the average 
background asymmetry (^A]^'^ are determined simultaneously as weighted averages over 
the accessible kinematic range. This method does not require an arbitrary selection of mass 
windows for the signal and background regions as in the classical side-band subtraction 
method. Moreover, it yields a smaller statistical error compared to the latter, reaching 
practically the lower bound of the unbinned likelihood method [25j. This is achieved by 
weighting every event with its analysing power aLL(Ar). The same procedure is applied 



dm dX 



a(j)n{s + b) 
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for A^J^ . The weighting factors are thus 

ws = P.f^^a^i. , = PJ-^D . (4) 
s + b s + b 

The target polarisation P^, as a time dependent quantity, is not included into the weights 
because including it may generate false asymmetries. Note that all events in the mass 
window -400 MeV/c^ < Mxn - M^o < 400 MeV/c^ of Fig. [1] are used. Since the factor 
s/{s + b) in ws vanishes for events far away from the central peak, these events do not 
contribute significantly to {Ag/g)^, but contribute to the determination of (^Aj^Y 

By considering sums over the different event samples eight equations are derived 
from Eq. ^ 

E^c. = «^(l + /?^/^) +7^(<)) (5) 



nt Pt,iWs,iWc,i t T,f* Pt,iWB,iWc,i 



for the two target cells before (t = u, d) and after (t = m', d') the target spin reversal, once 
weighted with ws and once with w-b (C = S, B). Here Nt is the number of events observed 
in cell t. These eight equations contain 10 unknowns which are (Ag/g)^, (^A'^'^ and eight 
acceptance factors = / a^4>^n^{s + b)wc dX. 

Assuming that possible acceptance variations affect the upstream and downstream 
cells in the same way, i.e. a^/aQ = aQ/afj, reduces the number of unknowns to eight. 
With an extra, much weaker assumption that signal and background events from the 
same target cell are affected in the same way by the acceptance variations, one arrives 
at a system of eight equations with seven unknowns. Possible deviations from the above 
assumptions may generate false asymmetries which are included in the systematic error. 
Using the set of eight equations (see Eq. (^), the gluon polarisation (Ag/g)^ and the 

background asymmetry (^A"^^ are determined with a standard least square minimisation 
procedure taking into account the statistical correlation between the number of events in 
a given target cell weighted by ws and by wb- The analysis is performed independently 
for the D* and D'^ samples. 

The quantities Pt, Pfj,, Qll and S/{S + B) are obtained as follows. For Pt, values 
averaged over about one hour of data taking are used, a timescale over which the assump- 
tion of a stable target polarisation was shown to be justified. The beam polarisation P^ 
is parameterised as a function of the momentum which is measured for each incoming 
muon. The photon-gluon analysing power, au[^{X)/D, is parameterised in terms of mea- 
sured kinematic variables. It depends on partonic variables not accessible experimentally 
and is obtained using a neural network [28j trained on a Monte Carlo sample for D* 
mesons. For this purpose PGF events were generated with AROMA [27] in leading order 
QCD, processed by GEANT to simulate the response of the detector and finally recon- 
structed like real events. It was checked that the MC simulation describes the background 
subtracted data distributions in z and sufficiently well. The scale, /x, used in the MC 
was chosen as the transverse mass of the produced charmed quark pair, and is sufficiently 
large to justify the perturbative approach. The correlation between the generated oll and 
the parameterised oll is 81% (see Fig. [2]). The same parameterisation is valid for the 
and the D* samples. 



6 



Finally, the signal purity, s/{s+b), as a function of the invariant mass for each event, 
is determined from a fit of the invariant mass distributions of the D* and samples. In 
this fit the signal is described by a Gaussian distribution. In the D* case the background 
function is the sum of an exponential and a Gaussian, the latter added to describe the 
refiection of the K tttt^ decay, where the 7r° meson is not observed. In the 

case the background is described by the sum of two exponential distributions. Note that 
not only the variation of the signal purity (or s/h) with the mass, but also with other 
characteristics of the event, is taken into account. This is achieved by a method |29] based 
on a multivariate approach starting with a parameterisation of the signal-to-background 
ratio, integrated over a window around the mass, {S/B)^^^. The window is of ±40 
MeV/c^ for the D* sample and ±30 MeV/c^ for the sample. The parametrisation is 
the product of 10 functions, each one depending on one of the 10 variables describing 
the event kinematics and the RICH response. Typically six bins are defined in each of 
the variables and the mass spectra are fitted in each bin of each variable to provide the 
values of the S/ B ratios using the signal and background functions described above. Each 
of the 10 variables is considered successively and the parameters of the corresponding 
function are adjusted to reproduce the S/B ratios in all bins in this variable. Adjusting 
the parameters for one variable affects the agreement obtained for previous variables and 
thus the adjustment procedure has to be repeated until convergence is reached and all 
S/ B ratios are reproduced simultaneously. 

Using this parametrisation, each sample {D* and D^) is split into intervals of 
{S/B)pg^j. and the mass spectrum is fitted separately in each of them. As an illustration the 
invariant mass spectra obtained in the highest interval of {S/ B)^^^ are compared in Fig. [T] 
to those obtained for the full samples. The signal purity for each event is obtained from the 
fit to the mass spectrum in the interval of (S'/i?)par containing the event and this value is 
adjusted to the exact value of {S/B)p^^ for this event. To validate the procedure the fit in 
each (S'/i?)par interval is integrated over the window around the mass peak to obtain the 
S/ B value and compared with the average value obtained from the parametrisation. The 
consistency obtained guarantees that using the (S'/-B)par in the event weights does not 
introduce a bias. In addition, it is checked that weighting the wrong-charge background 
(i^^vr+vTgJo^ and charge conjugates) with the parametrised values of the mass averaged 
signal purity, [S/ {S + -B)]par , does not generate any artificial peak at M^-k = M^o. 

5 Results for the gluon polarisation 

A value for (Ag/g)^ is obtained for each of the 40 weeks of data taking separately 
for the D° and the D* sample. The results (A5(/5()^=-0.421 ± 0.424(stat) for the 1)° and 
(Ag/ g)^— 0.541 ± 0.343(stat) for the D* sample, are the weighted mean of these values. 
The resulting background asymmetries, (^A^'^ = 0.003 ± 0.004 for the sample and 

Al'^) = 0.062±0.042 for the D* sample, are consistent with zero. Assuming that Ag/g{x) 
is approximately linearly dependent on x in the range covered, (Ag/g)^ gives a measure- 
ment of Ag/g{{x)), where (x) is calculated using the signal weights. This assumption is 
supported by the results of the COMPASS QCD analysis [5]. 

The major contributions to the systematic uncertainty are listed in Table [TJ The 
contributions from P^, Pt and / are discussed with more detail in Ref. [5J. To study 
the infiuence of false asymmetries, the final samples from Fig. [1] were subdivided into 
two samples using criteria related to the experimental apparatus, e.g. kaons going to the 
upper or to the lower spectrometer parts. The resulting asymmetries were found to be 
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Table 1: Systematic error contributions to (Ag/g)^ for D^{D*) channels. 



source 




source 


^(^) 


False asymmetry 


0.05(0.05) 


Beam polarisation 


0.02 


S/{S + B) 


0.07(0.01) 


Target polarisation Pt 


0.02 


All 


0.05(0.03) 


Dilution factor / 


0.02 




Total error 


0.11(0.07) 





compatible within their statistical accuracy, thus no false asymmetries were observed. An 
upper limit of the contribution of time dependent acceptance effects to the systematic 
uncertainty was estimated from the dispersion of the values for (Ag/g)^ and (^A^'j for 
the 40 weeks of data taking. Assuming that possible detector instabilities are similar 
for background and signal events and applying the method used in Ref. [5j leads to a 
conservative limit of 0.05 for both decay channels. 

Varying the procedure to build the parameterisation of s/{s + b), and in particular 
the functional form of the background fit, results in an error on (Ag/g)^ of 0.07 and 0.01 
for the D° and the D* sample, respectively. As expected, the uncertainty on s/{s + b) is 
larger for the D° case, where the signal-to-background ratio is smaller. To estimate the 
influence of the simulation parameters, i.e. charmed quark mass (varied from 1.3 GeV/c^ 
to 1.6 GeV/c^), parton distribution functions and scales (varied by a factor of 8), MC 
samples with different parameter sets were generated and oll was recalculated, resulting 
in an uncertainty on (Ag/g)^ of 0.05 and 0.03 for the and the D* sample, respectively. 
Other contributions, like radiative corrections and event migration between target cells, 
were studied and found to be negligible. 

The final value is the weighted mean of the two values for the D* and the sample 
and amounts to 

= -0.49±0.27(stat) ±0.11(syst) (7) 

in the range of 0.06 < x < 0.22 with (x) ^ 0.11, and a scale (/i^) ^ 13 (GeV/c)^. The 
contributions to the systematic uncertainty for each sample are added in quadrature to 
obtain the total error, 0.11 and 0.07 for the and D* sample, respectively. The larger 
value is chosen as a conservative estimate of the final error in Eq. ([7]). 

In Fig. [3] the above result is compared to other measurements of (Ag/g)^ and to 
two parametrisations from the NLO QCD analysis of the world data on the polarised 
structure function gi{x, Q"^), performed by COMPASS [5J: with AG >0 (broken line) and 
with AG* < (dotted line). The present result is consistent with previous measurements 
favouring small values of {Ag/g)^. Note that is the scale for the analysis of the SMC 
[13] measurement and the QCD analysis [5]. The scale of the present result is given by the 
transverse mass of the charmed quarks /i^ = ^ 13 (GeV/c)^. The other experimental 
points in Fig. [3] are given at /i^ ~ 3 (GeV/c)^. 

6 Asymmetry determination 

The data described in sections 2 and 3 also allow for the determination of the 
virtual photon asymmetry for production, A'^^^^ = ^m^^m contrast to 
(Ag/g)^ this asymmetry is independent of the interpretation in LO QCD. The asymmetry 
averaged over the full kinematic range would be largely diluted because of the large 
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dispersion of an. The asymmetry A^^'^^ ^ is thus extracted in bins of the transverse 
momentum of the with respect to the virtual photon, p!^ , and the energy of the 
in the laboratory system, -Edo. The bins were chosen such that the variation of a^^/D 
within each bin is small compared to the variation over the whole sample. In principle 
^7Af^D0x g^jgQ depends on the inclusive variables y and Q^, but an additional binning is 
not necessary because the dependence is very weak. This is clearly seen in LO, where 
j^^N^D X _ (^d-^^ I /\g I g _ In a given bin in and Ej^o the factor (qllZ-D) is almost 
independent of y and Q^, and the same is true for l^gj g. 

The asymmetry ^t^^^^^x jg obtained in exactly the same way as {/S.g/g)^, except 
that the factor qll is replaced by D in the definition of the signal weight in Eq. (jl]), 
i.e. ws = P^fDs/{s + b). This provides A'^^^^ ^{{Pt )A^do)) under the assumption 
that the bins of p^ and Eijo are small enough. It was verified that this approximation 
and the independence on y and are well fulfilled for the cross-section evaluated in LO 
QCD. At higher orders, the variation of the cross-section are expected to be similar and 
thus the approximations to remain valid. 

Table [2] gives A'^^^^°^ averaged over the and D* sample in each {pf^°,E£io) 
bin, together with the average of several kinematic variables. All averages are calculated 
with the weight ws = P^fDs/{s + h). The muon-nucleon asymmetry ^mJv^m'-D°x 
obtained by multiplying A^^^^ ^ by D{{X)). Both asymmetries can be used in global 
NLO QCD fits to constrain the values of Ag{x). 

As a cross-check we have calculated (Ag/g)^ from A"'^^^ ^ in each bin by dividing 
the asymmetry by the corresponding ollZ-D. Combining all bins we got a result consistent 
with the result in Eq. ([7]), with an increase of 5% in the statistical error. The contributions 
to the systematic error listed in Table [2] contribute as well to the systematic error of 
the asymmetries, except for the contribution of qll- This leads to a relative systematic 
uncertainty of 20% for A'^'^^^ ^ which is 100% correlated between the bins. 

7 Conclusion 

We have studied D° meson production in 160 GeV polarised muon scattering off 
a polarised deuteron target. The decays into pairs of charged K and tt mesons were 
selected using analysing the invariant mass distributions of identified Kvr pairs. Only one 
meson was demanded in each event. 

The data provide an average value of the gluon polarisation in the nucleon, (Ag/g)^, 
under the assumption that photon-gluon fusion to a cc pair is the underlying partonic 
process for open charm production, which is equivalent to a LO QCD approach. The result 
is {Ag/g)^= —0.49 ± 0.27(stat) ± O.ll(syst) at an average gluon momentum fraction, 
(x) ~ 0.11 and at a scale fi^ ~ 13 (GeV/c)^. This result is compatible with our previous 
result from the analysis of high-px hadron pairs but it is much less model dependent. 

The present measurement of the gluon polarisation in the nucleon, together with 
other measurements of COMPASS and HERMES, all situated around x ~ 0.1, point 
towards a small gluon polarisation at that value of x. This is a hint for a small value of 
the first moment, AG, of the gluon helicity distribution, although this in principle does 
not exclude a large value. 

The longitudinal cross-section asymmetries A'^^^^ ^ were also extracted from our 
data and are presented in bins of the transverse momentum and the laboratory energy of 
the D°. They may be used to constrain the values of Ag[x) in future global NLO QCD 
analyses. 
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Figure 1: Invariant mass distributions of the i^vr pairs for the D* sample (upper plot) and the 
sample (lower plot). The non-shaded histograms (left scale) show the total event samples 
while the shaded ones (right scale) show the events in the highest bin of {S/B)ps.T- 
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Figure 2: Correlation between the generated analysing power af^^ and the analysing 
power parameterised by neural network a^^ . 
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Figure 3: Compilation of the (Ag/g)^ measurements from open charm and high-px hadron 
pair production by COMPASS [i5\, SMC ^Mj and HERMES [13] as a function of x. The hor- 
izontal bars mark the range in x for each measurement, the vertical ones give the statistical 
precision and the total errors (if available). The open charm measurement is at a scale of about 
13 (GeV/c)^, other measurements at 3 (GeV/c)^. The curves display two parameterisations from 
the COMPASS QCD analysis at NLO [5] , with AG > (broken line) and with AG < (dotted 
line). 
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bin limits 




(y) 


iQ') 


iP?) 


(Ed) 


Dm) 




p?(GeV/c) 








(GcV/c)2 


(GcV/c) 


(GeV) 






0-0.3 


0-30 


-1.34 ±0.85 


0.47 


0.50 


0.19 


24.8 


0.57 


0.37 


0-0.3 


30-50 


-0.27 ±0.52 


0.58 


0.75 


0.20 


39.2 


0.70 


0.48 


0-0.3 


> 50 


-0.07 ±0.66 


0.67 


1.06 


0.20 


60.0 


0.80 


0.61 


0.3-0.7 


0-30 


-0.85 ±0.51 


0.47 


0.47 


0.50 


25.1 


0.56 


0.26 


0.3-0.7 


30-50 


0.09 ±0.29 


0.58 


0.65 


0.51 


39.4 


0.71 


0.34 


0.3-0.7 


> 50 


-0.20 ±0.37 


0.67 


0.68 


0.50 


59.6 


0.80 


0.46 


0.7-1 


0-30 


-0.47 ±0.56 


0.48 


0.53 


0.85 


25.2 


0.58 


0.13 


0.7-1 


30-50 


-0.49 ±0.32 


0.58 


0.66 


0.85 


39.1 


0.70 


0.17 


0.7-1 


> 50 


1.23 ±0.43 


0.68 


0.73 


0.84 


59.4 


0.81 


0.26 


1-1.5 


0-30 


-0.87 ±0.48 


0.50 


0.49 


1.21 


25.7 


0.60 


0.01 


1-1.5 


30-50 


-0.24 ±0.25 


0.60 


0.62 


1.22 


39.5 


0.73 


0.00 


1-1.5 


> 50 


-0.18 ±0.34 


0.69 


0.77 


1.22 


59.3 


0.83 


0.04 


> 1.5 


0-30 


0.83 ±0.71 


0.52 


0.51 


1.77 


26.2 


0.63 


-0.13 


> 1.5 


30-50 


0.18 ±0.28 


0.61 


0.68 


1.87 


40.0 


0.74 


-0.20 


> 1.5 


> 50 


0.44 ± 0.33 


0.71 


0.86 


1.94 


59.9 


0.84 


-0.24 



Table 2: The asymmetries A'^'^^^ in bins of p^" E^o for the and D* sample combined, together with the averages of several 
kinematic variables. Only the statistical errors are given. The relative systematic uncertainty is 20% which is 100% correlated between 
the bins. 



